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r-| \ Abstract 

We investigate whether the lightest scalar mesons a and k have a large tetraquark com- 
^ , ponent qqqq^ as is strongly supported by many phenomenological studies. A search for 

00 \ possible light tetraquark states with J^^ = 0^"*" and / = 0, 2, 1/2, 3/2 on the lattice is 

^ I presented. We perform the two-flavor dynamical simulation with Chirally Improved quarks 

O ' and the quenched simulation with overlap quarks, finding qualitative agreement between 

If) ■ both results. The spectrum is determined using the generalized eigenvalue method with a 

^ , number of tetraquark interpolators at the source and the sink, and we omit the disconnected 

contractions. The time- dependence of the eigenvalues at finite temporal extent of the lattice 
is explored also analytically. In all the channels, we unavoidably find lowest scattering states 
Ti{k)Ti{—k) or K{k)TT{—k) with back-to-back momentum /c = 0, 2ti / L , .. . However, we find 
an additional light state in the / = and / = 1/2 channels, which may be interpreted as 
the observed resonances a and k with a sizable tetraquark component. In the exotic repul- 
sive channels 1 = 2 and I = 3/2, where no resonance is observed, we find no light state in 
addition to the scattering states. 



% 



1 Introduction 

The only well established hadron states so far are mesons qq and baryons qqq. No exotic 
states like tetraquark [gg][gg], pentaquark qqqqq, hybrid qqG or molecular {qq){qq), {qq){qqq) 
have been confirmed beyond doubt, although there are several serious candidates in the light 
and hidden charm sectors. Perhaps the most prominent tetraquark candidate is the Z"'"(4430) 
resonance, discovered by Belle [T]: it decays to 7r~^il)', so it must have a minimal quark content 
ducc, but it has not been confirmed by Babar [2]. 

It is still not established whether the lightest scalar mesons a, n, ao(980) and /o(980) are 
conventional qq states or they have important qqqq or glue Fock components. A sizable glue 
component of isoscalar a is supported by some phenomenological studies [3], but we will 
not explore this Fock component in this work. We will focus on the qqqq Fock component, 
which arises in the case of a tetraquark [qq] [qq] or in the case of a mesonic molecule {qq){qq)- 
The tetraquarks [q'q'][gg] are composed of a scalar diquark (3c, f) and anti-diquark (3c,f) in 
L = 0; they form a flavor nonet and are expected to be light [HE]. A mesonic molecule 
{qq){qq) is composed of two color-singlet mesons (vr, K) held together by pion exchange [6]. 
Both qqqq interpretations expect that the 1 = 1 state [ussd) is heavier than the 7 = 1/2 
state {udds) due to m^ > tti^, in agreement with experimental ordering ^^^(ggo) > Tn,^. On 
the other hand, the conventional ud and us states can hardly explain the observed mass 
ordering. Both qqqq interpretations also naturally explain the large observed coupling of 
ao(980) and /o(980) to KK, which is due to the additional valence pair ss. 

In this paper we use a lattice QCD simulation to address the question whether the lightest 
scalar mesons a (/ = 0) and n [I = 1/2) have a sizable qqqq component. The quantities 
studied in our present simulation do not distinguish between a tetraquark and a mesonic 
molecule. When we use a word "tetraquark" below, we have in mind both types of exotic 
qqqq states. 

The a resonance is now widely accepted since its pole with rria = 441 ""^g MeV and 
Fg. = 54:4:~^ MeV was determined in a model- independent way [7j. The k resonance pole 
with ttIk = 658 ± 13 MeV and F^ = 557 ± 24 MeV was determined in a similar manner [8]. 
Both resonances have been recently experimentally confirmed [9j, but they remain slightly 
controversial. 

In order to extract the information about tetraquark states, lattice QCD simulations 
evaluate correlation functions with tetraquark interpolators at the source and the sink. In 
addition to possible tetraquarks, also the scattering states P1P2 (vrvr for 7 = 0,2 and Kit for 
/ = 1/2, 3/2) unavoidably contribute to the correlation function and this presents the main 
obstacle in extracting the information about tetraquarks. The scattering states Pi{k)P2{—k) 
at total momentum p = have discrete energy levels 



Ep^P^ ^ Ep^{k) + Ep^{-k) , with Ep{k) = \J m], + F and k = ^n (1) 

in the non-interacting approximation. The energy level Ep-^p^ ^ rrip-^^ + mp^ is low and 
makes an important contribution to the correlation functions. In order to identify possible 
tetraquarks, one has to extract several energy levels En in each isospin channel and then 



consider various criteria that could distinguish between the one-particle (tetraquark) states 
and the two-particle (scattering) states. We do not consider the more challenging 1 = 1 
channel, since there are two towers i^+(A;)A'°(— /c) and Ti{k)rj{—k) of scattering states. 

The lattice simulations [TOl [TTl [T2l [T3] have not yet provided the final answer to whether 
the lightest scalar mesons are tetraquarks or conventional qq mesons. All previous tetraquark 
simulations were quenched and they ignored disconnected contractions (cf. Fig. [3]). All 
simulations (except for [HI [12]) consider only the / = channel and the simulations [121113] 
extract only the ground state. The strongest indication for a as a tetraquark was obtained 
for rrtTT — 180—300 MeV in [10] by considering the lowest three energy level^j and the volume- 
dependence of the spectral weights. This impressive result on a meson was obtained from 
a single correlator using the sequential empirical Bayes method [10] and needs confirmation 
using a different method, for example the variational method used here. The first study 
that used the variational method to extract the ground and the excited energy levels in 
/ = 0, 1/2 channels was presented in [TTj, but the first excited state was found much higher 
than Ep^(^) + Ep^{—^). The reason for that was attributed to the unfortunate choice 
of the interpolators that had the same color and Dirac structure, while they differed only 
in spatial structure. For this reason we take interpolators with different color and Dirac 
structures in the present analysis, which enables us to extract the state Pi{^)P2{—^). We 
note that there have been few lattice simulation of tetraquarks or mesonic molecules in the 
related hidden charm sector [T^ . 

In this paper we determine a spectrum of states with J^*^ = O"^"*", p = and / = 
0, 2, 1/2, 3/2 on the lattice using the variational method with a number of tetraquark 
sources and sinks. We also determine the couplings (0|(9i|n) between the interpolators Oi 
and the physical states \n). This is the first dynamical simulation intended to look for 
tetraquarks and we also perform the quenched simulation in order to see whether there 
are any qualitative differences between the two cases. Our djTiamical simulation has two 
flavors of Chirally Improved quarks, while our quenched simulation uses overlap quarks. In 
this pioneering study, we are interested in the "pure" tetraquark states with four valence 
quarks qqqq and we prevent qqqq -H- qq 'H- vac ^ glue mixing by neglecting the disconnected 
contractions in J = 0, 1/2 channels, as in all previous tetraquark studieqj. 

The accurate lattice spectrum En as a function of the lattice size L in principle allows 
determination of the resonance mass and widths [TTJ [18] . The resonance appears as a state 
in the spectrum in addition to the discrete tower of scattering states. At the values of m.^L, 
where the non-interacting energies of the resonance and the scattering state would cross, the 
energy levels experience the largest energy shifts from the non-interacting values [T7[ [T8] . 
These energy shifts in principle allow the determination of the resonance width. In practice, 
the accurate determination of the excited energy levels on the lattice is challenging and only 
the p meson width has been reliably determined from the ground energy level in this way 
So far only simulations of the toy models were able to extract the scattering states and 



^The / = energy levels from [10] and the present work are compared in the Section 3.3. 
^The disconnected contractions have been recently taken into account in the study of the ground scattering 
states with / = [l5] and / = 1/2 [l6]. 



the width of the resonance from the ground and the excited energy levels [20] . 

In the present work, our excited energy levels are not accurate enough to allow for the 
determination of the a and k widths. We concentrate on a simpler and more realistic ques- 
tion: is there any light state in addition to Pi(0)P2(0) and Pi(^)P2(— x) ^'^ ^^^ attractive 
channels / = or / = 1/2? Such an additional state could be related to cr or k with sizable 
tetraquark components. Our main result for the spectrum, shown in Figs. [1] and [2|, indeed 
shows an additional light state and we propose a possible interpretation that it is a tetraquark 
state. The extracted energy of the additional state as a function of m,r qualitatively agrees 
with mo-,K('^7r) from unitarized Chiral Perturbation Theory (ChPT) [5]. 

The 0" and k are expected to become bound states at our heavier m.,^, where decays 
cr — 7- TTTT and n — )■ Ktt are no longer allowed kinematically [5]. In the range of ttItt, where this 
might occur, we find candidates for a and k close to the threshold, in qualitative agreement 
with the prediction of unitarized ChPT [5] and a lattice study of a toy-model with bound 
and scattering states [2T] . 

We also determine the spectra in repulsive channels 1 = 2 and / = 3/2, where no light 
resonance is experimentally observed. Our main purpose here is to verify that there is 
no light state in addition to Pi(0)P2(0) and Pi(^)P2(-^). The results in Figs. [U and [2] 
demonstrate that indeed we do not find any additional state. 

We explore two methods to distinguish the one-particle and two-particle states. The 
first is based on the time- dependence of the correlation functions and the eigenvalues of the 
variational method, which are explored analytically as well. The second method is based 
on the volume dependence of the couplings {0\Oi\n), which is explored in the quenched 
simulation and compared to the theoretical expectations. 

Some of our initial exploratory results have been published in proceedings [22] . 

We present the methods to extract the spectrum and the (0|Cj|n) couplings in Sect. 2. 
The analytical expectations for the time-dependence of the correlators and the eigenvalues at 
finite temporal extent are also given in this section, while certain derivations are delegated 
to the Appendix A. Our numerical results and their interpretation are given in Sect. 3, 
and we end with conclusions. In Appendix B we show that omission of the disconnected 
contractions cannot lead to an unphysical intermediate state with different isospin. 

2 Simulation 

2.1 Correlation matrix and interpolators 

We computed the same tetraquark correlation functions Cij{t) in a two- flavor dynamical 
simulation and in a quenched simulation. The purpose was to see whether there are any 
significant qualitative differences between the two cases. The results, presented in the Section 
[3l show qualitative agreement between dynamical and quenched results. The details about 
both simulation are given in Section 12. 4[ 

The tetraquark correlation function creates a state qqqq with chosen / and J^'" = 0~^~^ 
at t = and annihilates it at some later time t, where the projection to the total momentum 
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Figure 1: Tlie resulting spectrum E^ for / = 0, 2 in the dynamical (left) and the quenched (right) 
simulations. Note that there are two states (n = 1 and 2) close to each other in the / = case. 
The lines present the energies of non-interacting 7r(fc)7r(— A;) with k = N 271 /L and A^ = 0,1,^2. 
The plotted data is given in Tables [3] and HI 
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Figure 2: The resulting spectrum En for / = 1/2, 3/2 in the dynamical and the quenched 
simulations. Note that there are two states (n = 1 and 2) close to each other in / = 1/2 case. The 
lines present the energies of non-interacting K{k)7r{—k) with k = N 27r/L and A^ = 0, 1, v2. The 
plotted data is given in Tables El and SI 

p = is made. The time- dependence is obtained by inserting a complete set of physical 
states \n) with given quantum numbers 



a,{t) = {0\O,{t)O]m0)^^^, = ^(0|O.(x,t)O](0,0)|0) ^-^ J^ZfZf e 



-E„ t 



(2) 
with coupling Z" = (0|(9j|n) and temporal extent of the lattice T. The correlation matrix 
is used to extract energy levels En and couplings Z" for the tetraquark system. 

We consider five interpolators in the case of / = 0, 1/2 and three interpolators in the 
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case of / = 2, 3/2. The interpolators differ only in Dirac and color structure, but they have 
the same spatial structure: all quark fields are evaluated at the same space-time point. 

For 1 = and 1/2, where a resonances may be expected, the first three interpolators 
C^i=i,2,3 are products of two color-singlet currents (with sum over spatial components fi for 
Ci=2,3) and the flavor structure is derived in Appendix B. The last two interpolators are of 
well-known diquark anti-diquark type |1] 

Oit2,3 = E 2(rfrfn)(nrfrf) + i(nrfn)(nrfn) + i(<d)(<ci) - (nrfn)(Jrfd) , 

11=1,2,3 

/i=l,2,3 q=u,d,s 
0'=li' = [sTid^]a [u^Tid\a . (3) 

Here T = 7or'l"7o and 

Ti = 75 , r^ = 7^ r^ = 7'^75 , 14 = C75 , Ts = C , (4) 

while [g'^FQ] denotes a (pseudo) scalar diquark [q'^TQ]a = ^abcllb^Qc — Qj^Qc]- The / = 1/2 
tetraquark interpolators above are constructed to transform as |/, 13) = 1 1/2, 1/2) under 
SU{2)f and like su flavor state under SU{?>)f- 

We use a smaller three-dimensional interpolator basis for 1 = 2 and 3/2, which are not of 
our prime interest since no resonances have been observed experimentally in these repulsive 
channels. All these interpolators are of current- current type 

Ol:l2,,= E ^dT^^^)^dT>tu) 

M=l,2,3 

^Ka= E (^"r^)(^"rf..) (5) 

M=l,2,3 

with Fj as defined in (JH). 

Fig. [3] shows the contractions which enter the correlation matrix with our tetraquark 
interpolators at the source and the sink: / = has all three contractions, 1 = 1/2 has 
contractions (a,b) and / = 2, 3/2 have only the connected contraction (a). In this pioneering 
study, we are interested in physical states with four valence quarks qqqq ( "pure" tetraquark 
states and P1P2 states) and we therefore take into account only the connected contractions. 
The singly (doubly) disconnected contraction couples also to qq (vacuum and glueball) states 
and we ignore them in order to be able to attribute a definite valence qqqq quark structure to 
the obtained physical states. Another reason for omitting the disconnected contractions is 
that they are difficult to evaluate and they are often noisy. We note that it is not legitimate to 
ignore the disconnected contractions in a proper field theory, as it leads to the violation of the 
unitarity. A possible effect of this approximation on our results is discussed in Section 1X51 and 
in Appendix B. We choose this approximation in our study with the given physical motivation 
and leave a proper study of physical states including mixing qqqq O gg O vac O glue for 
the future. 
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Figure 3: Quark contractions for our tetraquark correlators: connected (a), singly disconnected (b) 
and doubly disconnected (c). Only connected contractions are taken into account in our simulation, 
for reasons explained in the text. 



2.2 Variational method at T — > oo 

The extraction of the energies En and the couphngs Z^ from the correlation functions ([2]) 
using multi-exponential fits is unstable. Instead, we use our N x N correlation matrix to 
compute the eigenvalues A"'(t) and eigenvectors u"(t) of the generalized eigenvalues problem 
|23l|24] 

C{t)ir{t) = X^{t,to)C{to)u^it) . (6) 

The energy can be extracted from 



A"(t) *-^ e- 



E„{t-to) 



(7) 



when the temporal extent of the lattice T is very large. The error on the extracted energy 
En due to finite basis A^ is (^(e-^^^+i"-^")*) for to < ^ < 2to |21]- We will demonstrate that 
our results are almost independent of to for to ^ [I5 4], while they get noisier for to > 5. Our 
main analysis is based on t G [6, 10], so the condition t < 2to is satisfied or close to being 
satisfied. The exponential time-dependence ([7]) would apply only for infinite temporal extent 
T and we discuss the significant effect of finite T on Cij(t) and A"'(t) in Sect. 12.31 

The eigenvectors u^{t), which satisfy the orthogonality relation {'iP,C{t)u^"^^) oc 6nm, 
allow us to determine the couplings Z" at large t (see for example the derivation in 



\{Omn)\ 



E.a.(t)n^(t)| 



,E„t/2 



VE 



n* 

Im I"; 



(t) cut) u^{t)\ 



Note that the normalization of 'u"(t), which is arbitrary, cancels in ([8]). The error on the 
extracted coupling Zf due to the finite basis A^ is 0{e~'^^^+^~^"''^°) for t < 2 to and for fixed 
t — to [21]. Extracting the ratio of couplings for a given state \n) to two different interpolators 
is particularly straightforward 



zr 



Zn 



\J:kC^kit) Ul{t)\ 

Ek'CMtK'it)\ 



Extracting |Zf | itself 



E.c'..(t)<(t)| 



VTZ\yrWcut)KM 



\z?\ 



-Bt 



(9) 



(10) 



requires fitting the LHS to the form on the RHS. We always verify that the fitted coefficient 
B in the exponent is consistent with -E'„,/2 as obtained from A"'(t). 



2.3 Effect of finite T on tiie correlation matrix and the eigenvalues 

Our temporal extents (T = 32 for the dynamical and T = 28 for the quenched simulation) are 
not very large, so we need to understand the effect of finite T on Cij{t) and A"(t). At finite 
temporal extent T, the time-dependence e~^"* gets modified depending on the boundary 
conditions and the nature of the states. We use anti-periodic boundary conditions in the 
time direction for quarks and anti-quarks. 

If a single one-particle state \n) dominates the correlator, the diagonal correlator Cuit) 
behaves as 

Q^^t) '^^' |Zr|2 [e-^-* + e-^"(^-*)] . (11) 

If the correlator is dominated by a two-particle state \n) = \P1P2), it behaves as (see Ap- 
pendix A of this paper and Appendix A of [11] and [26] ) 

C,,{t) '"^^ iZrP [e-^-* + e-^"(^-*)] + |Zr|2 [e-^'p.te-^^^iT-t) + g-mp,tg-mp,(T-t)] ^ ^^3) 

where En is two-particle energy. Let us consider the relative importance of the couplings 
Zf = (OlCilPiPa) and If = (P/|Ci|P2) [H], which will be needed in our further study. Both 
matrix elements have similar structure, therefore one expects that Z^ is of the same order of 
magnitude as Zf . However we do noi expect that (Pj' |Cj|P2) is exactly equal to (0|Oj|PiP2): 
in the second case Oi annihilates the (interacting) state P1P2, where Pi and P2 existed at 
the same time and therefore interacted; in the first case Oi annihilates the P2 and creates 
P/, so Pi and P2 never exist at the same time and there is no interaction between them. We 
therefore believe that in the interacting theory (P/|(9j|P2) 7^ {0\Oi\PiP2) and Z^ ^ Z^, but 
we expect they are of the same order of magnitude. 

In reality, several physical states contribute to the correlation matrix and the time- 
dependence of the eigenvalues becomes more complicated. We are in particular interested 
in the cases where two-particle states and also possible one-particle (tetraquark) states con- 
tribute to the correlation matrix. We study the generalized eigenvalue problem for such a 
situation in the Appendix A. We find that 

(i) the eigenvalue corresponding to the one-particle state would have a time-dependence 

proportional to e~-^"* -|- e"-^"^^"*) and 

(ii) the eigenvalue corresponding to the two-particle state would have a time-dependence 

proportional to e"^"* + e-^"(^-*) + /?[e-'^^i*e-'"^2('^-*) + e-'"^2*e-'"^i(^-*)] with two- 
particle energy En 

only if the following relations would apply exactly: (P/|Oj|P2) = (0|Oi|PiP2) or more gen- 
erally if (Pi |0j|P2) = R {0\Oi\PiP2) with R independent of i. However, we argued in the 
previous paragraph that these relations do not apply exactly and so the time-dependence 



of eigenvalues is more complicated than in statements (i,ii) above. In the Appendix A we 
argue (although this has not been rigorously proved) that the two-particle as well as the 
one-particle eigenvalues have the general formo 

A"(t) = W" [e-^-* + e-^"(^-*)] + W" [e-™,p^te-mP2(T-t) ^ ^-rap^t^-mp^{T-t)^ ^ ^^3^ 

In our analysis, we extract En from eigenvalues A"(t) using a three-parameter fit {En-, w, w), 
where masses of Pi 2 = n, K are fixeqj to the measured values given in Tables [1] and [2j In 
case of equality {P{\Oi\P2) = {0\Oi\PiP2), we would expect w = ior one-particle state and 
w = w for two-particle state. We note that the extracted w" and w" in Tables [3] and H] are 
of similar magnitude. 

Finally we note that the form (fT3|l generally applies only in case when all diagonal and 
non-diagonal correlators are symmetric with respect to t — > T — t and that all our Cij{t) have 
this property. In case when some parts of non-diagonal correlators are anti-symmetric with 
respect to t ^ T — t, "backward propagation" described in Section IIF of [27] may appear. 



2.4 Details of dynamical and quenched lattice simulations 

Our correlation functions are constructed based on gauge configurations and quark propa- 
gators for two cases: the two-flavor dynamical simulation [28] and the quenched simulation 

m- 

• Details of the two-flavor dynamical simulation together with the ground state 
hadron spectroscopy are given in [28], while some related results based on a similar 
setup are given in [271 129]. In that simulation the Liischer-Weisz gauge action [30] 
and two- flavors of dynamical degenerate Chirally Improved (CI) quarks [31] have been 
used. The gauge fields are periodic in all four space-time directions, the fermion field 
anti-periodic in the time direction. One level of stout smearing [52] is applied to the 
gauge configurations, which is considered as a part of full Dirac operator. We use three 
ensembles (C,B,A) with the same lattice volume 16^ x 32 and three different u/d quark 
masses, corresponding to m,r — 318 — 526 MeV (see Tabled]). The lattice spacings, also 
given in Tabled! have been determined using rg = 0.48 fm and are close to a ~ 0.15 
fm for all three ensembles. 

The valence u, d, s quarks are also of Chirally Improved type. The valence u/d quark 
masses are always fixed to the dynamical u/d quark masses in our study. The valence 
strange quark masses is fixed from rrin. 

All quark sources and sinks are Jacobi-smeared by applying 

[UMt)6,^ly + U^^{x-j,t)6,_jJ , (14) 

n=0 j=l 

■^This form applies if only one two-particle state P1P2 makes significant contribution at t ~ T/2. 
^We verified that the variation of the results is negligible if m^^/^ are varied in the ranges given in Tables 
[T]and[l 
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Figure 4: Cosh- type effective mass (fT5]) for the diagonal correlators Cii(t) at rriTr 
the dynamical simulation. The lines indicate 2771,^ or m.,^ + rnx- 



469 MeV in 



which is invariant under rotations and preserves the interpolator quantum numbers. 
For this analysis we use a single ( "narrow" ) smearing with the values of A^ and k in 
Table [H which are chosen to give the source/sink Gaussian width of approximately 
0.27 fm. 



Details of the quenched simulation are presented in [TOl [33]. It employs overlap 
valence quarks, which have exact chiral symmetry even at finite lattice spacing. The 
gauge fields are generated using Iwasaki actions, where the lattice spacing a = 0.200(3) 
fm is determined using J't^. Our main analysis is based on the volume 16'^ x 28, while 
a smaller volume 12^ x 28 with the same lattice spacing is used for the study of the 
volume-dependence of couplings {0\Oi\n). All u, d, s quarks have point-like sources and 
sinks, while the strange quark mass is fixed from m^. We use three u/d quark masses, 
corresponding to m-„ = 230 — 478 MeV (see Table [2]). 
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Figure 5: The typical cosli-type effective masses according to (fTSj) (open symbols) and the cor- 
rected effective masses (|16p (filled symbols) for two lightest states in / = 0, 1/2 channels and for 
the lightest state in / = 2, 3/2 channels. The full 5x5 correlation matrix is used for / = 0, 1/2 
and the 3x3 correlation matrix is used for 1 = 2, 3/2. The figure shows as an example our results 
for the dynamical simulation with m-,^ = 469 MeV and to = 3, with w fitted via ()13p in the range 
t G [8, 15]. The lines represent 2771,^ (/ = 0, 2) and m^ + rriK {I = 1/2, 3/2). 



3 Results 

3.1 Time-dependence of diagonal correlators 

The effective masses for diagonal correlators Cuit) for four isospins are displayed in Fig. HI 
We show the cosh-type effective mass, defined as 



m 



<^ff + e 



e// 



(T-t) 



Fit 



e "ff 



(t+i) 



, -m' , (T-t-1) 



F{t) = Cu{t) or r{t) 



(15) 



so that JTig// equals the energy E ii F{t) = w[e~^^ + e~^^'^~''^^. The observed effective masses 
have sizable excited state contributions at small t and they have significant drop for t > 10, 
which indicates that a two-particle state \n) dominates Cii{t) at t ~ T/2 for each i and / 
(see Sect. [23]). 

3.2 Time-dependence of eigenvalues 

Due to the lack of reliable plateaus in the diagonal elements of the correlation function 
we turn to use the eigenvalues of the variational method in order to extract the spectrum. 
Typical cosh- type effective masses flTSjl for eigenvalues are plotted by empty symbols in Fig. |5l 
The excited state contribution at small t is smaller than for the diagonal correlators and rrieff 
reach a short plateau, but then they drop significantly for t > 10. We note that the drop is 
less significant for smaller 771^^ and more significant for larger m.,r, as illustrated in [11]. The 
drop demonstrates that eigenvalues don't have a simple e~^"^ + e"-^"*^-^"*) time-dependence. 
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Figure 6: Typical effective masses of the eigenvalues A"(t) for 1 = 0, 2. Corresponding ratios 
\Z^/Z^\ at given n are also shown (j is the largest component). The full 5x5 matrix is used for 
/ = and the 3x3 correlation matrix is used for 1 = 2. Results for a specific u/d masses and 
to = 1 in the dynamical simulation are shown. The lines present the energies of non- interacting 
7r(A;)7r(-A;) with k = NItt/L and A^ = 0, 1, ^2. 

but have a more complicated time-dependence flT3l) due to the presence of two-particle states 
TTTT or Kir in the box with finite T. The corrected effective mass takes that effect into account 



F{t) - w[e-"'Pi'e-"'P2(T-t) ^{t^T-t}] 
F{t+1) -^Z'[e-™^i(*+^)e-'^^2(^-*-^) + {t^T-t}] 



e "''^ff^ + e" 



'<!// 



(T-t) 



'''eff 



(t+l) 






(T-t-1) 



(16) 



so nieff is equal to En when F{t) = Cuit) or F{t) = A"(t) has the form ( fT3l) . Such me// is 
obtained after w has been determined by fitting A"(t) to ( TT3|) . The corrected effective mass 
is presented by the full symbols in Fig. [5l it is flat and it demonstrates that eigenvalues really 
have time-dependence of the form (TT3|) . 

We will always use the fitting form (fT3|) and the corrected effective mass ( 1X6]) for A"~q'-^ ,2 

and for A"^2 3/2- -^^^ higher states A^^qj^^ ^^^ -^"=23/2' ^^^ error-bars on A"(t) are large 
at t > 10, where the finite T effect is significant, such that the three-parameter fit with 
{En, w"-, w^) is not stable. In this case we will present cosh-type effective masses fll5j) and 
we will fit to 

X^{t) = w" [e"^"* + e-^"(^-*)] (17) 



at rather small t, where the finite T effect is not significant. 

Let us note that the significant effect at finite T = 28, 32 prevents us from a reliable 
determination of the energy shifts AEn = En — mp^ — mp^ , which require very long stable 
plateaus. Therefore, we do not aim at determining the energy shifts, but we determine the 
spectrum -E'„=i,2,3 itself with a reasonable precision. 
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3.3 Energy levels £■„ and couplings Z" 

/ = and 1 = 2 

The typical effective masses for J = 0, 2 are collected in Fig. [61 The ratios of couplings 
\Z^/Z^\ extracted via (|9]) are also shown. The lines display the three lowest energies of 
'K{k)'K{—k) in the non-interacting case. 

In J = case, we find one state with energy close to 7r(0)7r(0), another state with energy 
close to '^(^)'^{—'^) and we also find an additional light state (close to the lowest state). 
This applies for all quark masses and for the dynamical as well as the quenched simulation. 
We varied to ^ [1;4] and performed the diagonalization of 5 x 5 and all possible 4x4 and 
3x3 sub-matrices. We find that the extracted energiecl E^ and couplingqj |Z"| are almost 
independent of these choices for all quark masses and both simulations, as demonstrated for 
a specific case in Fig. [3 

In the 1 = 2 case, we find one state with energy close to 7r(0)7r(0), another state with 
energy close to 7r(^)7r(— ^) and no additional light state (see Fig. |6]). Again, this applies for 
all quark masses, both simulations and for the range of to ^ [I5 4]. In this case we use only a 
3x3 matrix ([5]), which is probably not large enough to capture the energy of 7r(^)7r(— ^) 
exactly (it naturally comes out too high). We point out that our intention was not to capture 
the energy of 7r(^)7r(— ^) correctly, but to verify that there is no light state in addition to 
7r(0)7r(0) in the 1 = 2 channel. 

The final result for the dependence of the extracted spectrum on m^r in both simulations 
is presented in Fig. [TJ Tables E] and H] provide the corresponding numerical results, together 
with the choices of to, interpolator sets and fit ranges. The values are obtained using an 
uncorrelated fit to f lT3|) or flTTj) and the error-bars are obtained with the single elimination 
jack-knife method. 

Finally, we compare our / = results to those obtained from a single tttt correlator using 
the sequential empirical Bayes method (SEBM) [10]. Above 300 MeV, we get degenerate 
results for the two states below 7r(^)7r(— ^). For m^ > 300 MeV, the authors of [lO] were 
unable to separate these two states using the SEBM; there is no disagreement since the 
SEBM is not designed to resolve degenerate states. Although the rriT^ > 300 MeV results 
were not included in the publication, the authors of [TO] do present results for pion masses 
in the range 182 MeV to 250 MeV. Both groups have analyzed a pion mass of 230 MeV. 
For this, the authors of [TU] get the first state close to 7r(0)7r(0) and the third state close 
to 7r(^)7r(— ^), just as we do. Although the detailed results are somewhat different, both 
groups obtain a second state below 7r(^)7r(— ^), which is the crucial result. 



^As an exception, all energies in Figs. [7] and |9] are extracted using the fit form ((T7)) at rather small 
t G [7,10], where a finite T effect is not significant. The resulting En are just intended to demonstrate 
independence on to and on the choice of the interpolator set. 

^The |Z"| was determined via (fTO]) by fitting in the time range indicated in the plot. 



13 



1=0, dyn. simulation, m =469 MeV, fit t=[7,10] 
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Figure 7: Extracted energies En and couplings Z" = (OlO^ln) in the / = channel for various 
choices of to and interpolator sets. A typical case with m-jr = 469 MeV in the dynamical simulation 
is shown. On the horizontal axis choice = 1,3,5,7 correspond to sets 012345, O12451 Oi345> 0123, 
respectively, where to = 2 is taken in all these cases. Similarly, choice = 2, 4, 6, 8 correspond to the 
same four sets and to = 3. The values E3 from interpolator set 01.2,3 are naturally too high since 
a 3 X 3 matrix cannot reliably provide E3. The Z^~ ' for the ground and the first excited states 
are independent of the choice to a good precision. For the second excited state, we show only the 
largest component Z2^^, which is expected to be the most reliable among all components. 



/ = 1/2 and / = 3/2 



The effective masses and ratios |Zf/Z"| for I = 1/2, 3/2 are shown in Fig. |8l while the 
resulting spectrum E^ is given in Fig. [2] and in Tables [3l HI 

The conclusions regarding I = 1/2 are completely analogous to the / = case above: 
there is a light state in addition to ^(0)71(0) and K{^)7r{—^). Results for the exotic 
/ = 3/2 channel are analogous to results for 1 = 2: there is no light state in addition to 
J^(0)vr(0) and i^(^)vr(-f ). 

This applies to all quark masses, to both simulations and to any choice of to G [1, 4]. We 
performed the diagonalization of the 5x5 and all possible 4x4 and 3x3 sub-matrices 
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Figure 8: Analogous to Fig. [H but for / = 1/2 and / = 3/2. The lines present the energies of 
non-interacting K{k)iT{ — k) with k = N2'k/L and A^ = 0, 1, \/2- 



1=1/2 , dyn. simulation, m =469 MeV, fit t=[7,10] 
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Figure 9: Analogous to Fig. [71 but for I = 1/2. 
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Figure 10: Comparison of Zf{L = 16) with Zf{L = 12) and (l|)3/2zf (L = 12) for the ground 
state n = 1, I = 0, 2, 3/2 and ttt-tt = 342 MeV in the quenched simulation. The couplings for the 
excited states n>2 have sizable errors and are not shown. All couplings are obtained via (jlOp by 
a fit in the range t E [7, 10] from a 3 x 3 matrix (/ = 2, 3/2) or a 5 x 5 matrix (/ = 0) at tg = 1- 



in / = 1/2 case. We find that the extracted E^ and |Z"| are almost independent of these 
choices for all quark masses and both simulations, as demonstrated for a specific case in 
Fig.l 



3.4 Volume dependence of Z, 



For completeness we provide now the volume dependence of couplings Z^ in the case of the 
quenched simulation, which was performed at two volumes 16'^ x 28 and 12^ x 28 at the same 
lattice spacing a = 0.200(3) fm. We are unable to show the analogous volume dependence 
in the case of dynamical simulation as it was performed on a single volume. 
The expectation for the L— dependence of Zf{L) is [101 [HI [IE] 

• Zf (16) ~ (i|)3/2^f (12) ~ 0.65 Zr(12) in case when \n) is two-particle state P1P2 

• Z"(16) ~ Z"(12) in case when \n) is a one-particle state (resonance) 

but these two behaviors are observed in practice only when eigenstates have very long and 
stable plateaus, as pointed out in [Mj . 

Figures [TO] and [n] compare Zf(16) with Zf(12) and (i|)3/2Zf (12) for all J. We have 
verified that the dependence of Z" on the choice of to and interpolator set is well below the 
(sizable) error-bars on Z^ for all i, n, I. 

The ground state n = 1 for 1 = 2, 3/2 is expected to be tttt or Kn since there are no 
observed light resonances in these repulsive channels. We determined couplings Z^^^ from 
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Figure 11: Comparison of Zf (L = 16) with Zl'{L = 12) and {^f/'^Z^{L = 12) for the three 
states n = 1,2,3 and / = 1/2 at 111^^ = 342 MeV in the quenched simulation. For the second excited 
state n = 3, we show only the largest component Z2~^, which is expected to be the most reliable. 
Couplings Z are obtained via pU|) by fit in the range t G [7, 10] from 5x5 matrix Q and to = 1. 



the variational method flTOj) as well as from the diagonal correlators Cii{t) using a fit flT2l) 
and both results agree within the error-bars. Our observed couplings seem to be roughly 
consistent with Z^{16) ^ (y|)^^^^[''(12), which applies for two-particle states. The relation 
is however not satisfied exactly, which is not surprising given our short plateaus and sizable 
finite T effect, as pointed out in [51] . 

The error-bars on extracted Z^ for / = 0, 1/2 in Figs. [TD] and [TT] are to large to distinguish 
between the one- and two-particle stateCl- 



3.5 Interpretation of the results 

We now discuss our interpretation of the observed spectrum in Figs. [1] and O which has the 
same distinctive features in the quenched and dynamical simulations. 

1 = 2 and / = 3/2 

Let us first focus on the repulsive channels 1 = 2 and / = 3/2, where n = 1 and n = 2 states 
are well separated in energy. We interpret the n = 1 ground state as 7r(0)7r(0) (i^(0)7r(0)) 
since it has its energy close to 2m^ (??7,^ + rrix) for / = 2 (/ = 3/2). The time dependence 
of the diagonal correlators Cii{t) in Fig. HJcan also distinguish the one-particle fllip or two- 



^The errors on Z" ' for / = 0, 1/2 in the dynamical simulation are smaller than in the quenched 
simulation (see Figs. [7] and [9]). 
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particle (TT2|) nature of the state n = 1, that dominates the correlatocl. The drop of cosh-type 
effective mass near t ^ T/2 in Fig. H] speaks in favor of two-particle state n = 1. The volume 
dependence of Z^^^ couplings in Fig. ^U\ is also roughly consistent with expectation for a 
two-particle state. The most important feature is a large gap between the n = 1 and n = 2 
states, so we do not observe any light resonance in / = 2, 3/2 channel, as expected for these 
repulsive channels. We interpret the n = 2 state as Pi(^)P2(~x) scattering state. It has 
an energy somewhat above Ep^{^) + Ep^{—^) ([T]), which is most probably due to the small 
3x3 interpolator basis, which is not big enough to capture -E'„=2 energy well. 

/ = and / = 1/2 

Now we turn to the more interesting attractive channels 1 = and 1/2, where broad reso- 
nances a and k may be expected. In each channel we observe two light states (n = 1,2) near 
the threshold mp-^ + mp^ and a third state n = 3 nicely consistent with Ep^(^) + Ep^{—^) 
(OP). It is natural to interpret the n = 3 state as the Pi(^)P2(— x) state, while the most 
interesting question is the nature of n = 1,2 states. 

States n = 1 and n = 2 are orthogonal to each other according to {u^,C{t)u^ ) oc (5„„/, 
so they must correspond to two distinct physical states. There is only one scattering state 
with given |/, J3) = |0, 0) in this energy range, and there is no way to envisage that the other 
state may also be a scattering state. This leads to a possible interpretation that the other 
state corresponds to a resonance o" for J = and resonance k for / = 1/2. 

Before considering this interpretation in more detail, let us point out again our most 
severe approximatioiij, which amounts to omitting the disconnected contractions in Fig. |3l 
One might wonder if this omission could lead to the unphysical appearance of an additional 
scattering state with wrong flavor near threshold. In the Appendix B we explicitly show that 
the scattering state |/ = 2,1^ = 0) cannot enter as an intermediate state in our connected 
correlator {O^^^IO"^^^^) , since the connected part of the matrix element (2, 0|(9^^^°) vanishes. 
The |/ = 1, Is = 0) state could also not appear as an intermediate state with E ~ 2mT^, since 
there is no 7r(0)7r(0) state with J^ = 0~^ and |/ = l,/3 = 0). In Appendix B we also show 
that the |/ = |,/3 = |) state cannot enter as intermediate states in our connected / = 1/2 
correlator, since the connected part of the matrix element (|, ^\0'^^^^^'^) vanishes. Therefore 
we believe that the omission of the disconnected contractions cannot be responsible for the 
appearance of the additional light scattering state. Our results with two light states in the 
/ = 0, 1/2 channels stimulate a future lattice simulation to search for low-lying states in 
1 = 0, 1/2 channels using the connected as well as the disconnected contractions. 

After these cautionary remarks, we examine the interpretation that one of the low-lying 
states is a scattering state and the other is a for 7 = and k, for / = 1/2. According to this 
interpretation, the resonances a/n found in the simulation are pure qqqq states and have 



^Thc n = 2 state is much heavier in case of / = 2, 3/2 and it dies out near t ~ T/2. 

^The other approximations are finite volume and lattice spacing, absence of dynamical strange quark 
and unphysical u/d quark masses, but these approximations could not lead to the presence of unphysical 
scattering state near threshold. 
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no qq Fock component: we are able to attribute definite Fock component to the simulated 
states since we used tetraquark interpolators and omitted the singly and doubly disconnected 
contractions in Fig. El The physical resonances (t/k, that correspond to the simulated states 
in Nature, therefore have a qqqq Fock component (probed in our simulation), but they 
may also have an additional qq component (not probed in our simulation without singly 
disconnected contractions). 

We already mentioned that the presence of the additional light state has to be confirmed 
in a future simulation, which takes into account also the disconnected contractions. Such 
a simulation will however not be able to claim the presence of the qqqq Fock component 
due to qqqq -^ qq mixing via singly disconnected contractions. Once the additional light 
state is confirmed in a simulation with disconnected contractions, the results of our present 
simulation will make case for the presence of qqqq Fock components in a and k. 

Now, let us attempt to establish which one of n = 1, 2 states is a candidate for the one- 
particle state a or k. Based on the drop of the cosh- type effective mass for Cait) near t ~ T/2 
(Fig. H]) we expect that a state \n) that dominates Cait) has a two-particle nature (IT^ . Our 
states ra = 1,2 are close to degenerate, so Cait) seems to be dominated by the state n = 1 
that has bigger coupling |Z""^| > \Z^^'^\ (Figs. \7\ [9l [TT]). This would indicate that ra = 1 is 
the scattering state Pi(0)P2(0), while n = 2 corresponds to a or k, although we cannot claim 
that with complete certainty. The errors on the couplings Z'^{L = 12, 16) for / = 0, 1/2 are 
too large to distinguish one- and two- particle behavior based on L— dependence of Z", as 
noted in Section 13.41 

Our interpretation that E""^^ corresponds to the scattering state Pi(0)P2(0) and E^^^ = 
mcr^K corresponds to a and k is in agreement with several expectations: 

• The a and k are resonances for our lower niT^ and so they are expected to be above 
the Pi(0)P2(0) scattering states, which is supported by the energies in Figs. [Hand [2] 

at low ?7i^. 

• The a and k are expected to be bound states for our higher ?7i^ [5]. In the range of m.^ 
where this might occur, we find a and k, almost degenerate with Pi(0)P2(0), which is 
in qualitative agreement with expectation from the lattice study of a toy model with 
loosely bound and scattering states [2T]. The authors [21] show that a loosely bound 
state lies slightly below the rrip-^ + nip^ and the Pi(0)P2(0) scattering state slightly 
abovqlj mp^ + rap^, but our study is not accurate enough to reliably extract these 
small energy shifts. 

• Our resulting m^r- dependence of rrifj^^ijnT,) = E'^^^{m.„) is in qualitative agreement 
with the prediction of ?7i(j,K('7i7r) within unitarized ChPT [5]. The authors [5] predict 
that ??7,o-,K are very close to threshold for mT^/m^y G [2,3], which is supported by our 
results in Figs. [1] and [2j The authors [5] find that a and k transform from a resonance 
to a bound state at about m^ ~ 350 MeV, which corresponds to intermediate 771,^ in 



^°The lowest scattering state Pi (0)^2(0) in the attractive channel is expected to be below rap^ + mp^ if 
there is no bound state. 
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our simulation. We note that at m.,^ > m^y used in our simulation, the a and k in 
are expected to be significantly narrower than the broad resonances observed in the 
experiment [S]. 

4 Conclusions and outlook 

We determined the energy spectrum and the couplings (0|(9i|n) for the states with J^'-' = 
O"*"^, isospin / = 0, 1/2, 3/2, 2 and p = using a number of tetraquark interpolators qqqq 
at the source and the sink. We omitted the disconnected contractions in Fig. [31 Our main 
question is whether there are any light states in addition to the towers of scattering states 
Pi{k)P2{—k) with k = 0,27r/L, .. and P1P2 = vrvr or Kir. Such additional states could be 
related to resonances a or k with a sizable tetraquark component. 

The resulting spectra in Figs. [1] and |2] show qualitative agreement between the dynamical 
and the quenched simulations. In the repulsive channels 1 = 2, 3/2, where no resonance 
is expected, we indeed find only the scattering states Pi(0)P2(0) and Pi(^)P2(— x) '^i^h 
no additional light state. In the attractive channels 1 = 0, 1/2 we find two (orthogonal) 
states close to the threshold mp-^ + rup^ and another state consistent with Pi(^)P2(— ^), 
so we do find an additional light state. This leads to a possible interpretation that one 
of the two light states is a scattering state Pi(0)P2(0) and the other one corresponds to a 
resonance a for / = and resonance k for / = 1/2. According to this interpretation, the 
physical resonances a and k have a nonzero qqqq Fock component, since the corresponding 
states in our simulation couple to the tetraquark interpolators. Along these lines, these 
physical resonances could not be pure qq resonances since pure qq resonances do not couple 
to tetraquark interpolators in absence of singly disconnected diagrams in our simulation. Our 
candidates for a and k have an rrij^ dependence in qualitative agreement with expectation 
from unitarized ChPT [5]. 

The volume dependence of the couplings (0|(!?i|n) for 1 = 2, 3/2 is roughly consistent with 
the above interpretation, while errors on J = 0, 1/2 couplings are too large to distinguish the 
one- and two-particle states based on this criterion. We also use the time- dependence of the 
correlators and eigenvalues at finite temporal extent as criteria for distinguishing the one or 
two particle states, and we demonstrate that one of the two light states in the / = 0, 1/2 
channels is a scattering state. Along the way, we derive the analytical conditions, which have 
to be satisfied so that the eigenvalues corresponding to the one-particle states would have a 
simple time dependence proportional to e~^"* -|- e~^"^'^~^\ 

We explored the possibility whether the omission of the disconnected contractions could 
lead to "unphysical" light eigenstate with J = 2 in the 1 = channel. We explicitly verified 
that 1 = 2 state cannot enter as intermediate state in our connected 1 = correlator. Simi- 
larly, I = 3/2 state cannot enter as intermediate state in our connected J = 1/2 correlator. 

The ultimate method to study a and n on the lattice would involve the study of the 
spectrum and couplings in presence of the disconnected contractions and the qqqq ^ qq -^ 
vac f-)- glue mixing, using interpolators that cover these Fock components. The recently 
proposed distillation method [35] could prove useful for determining the correlators with 
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Pi{k)P2{—k) interpolators or disconnected contractions. Such a study has to be done as a 
function of lattice size L in order to extract the resonance mass and width using the Liischer's 
finite volume method (TTJ dSl 1^ . 
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Appendix A: Effect of finite T on eigenvalues of the 
generalized eigenvalue method 

Here we consider the effect of finite T on the eigenvalues of the generalized eigenvalue prob- 
lem. We take a simple example with only two physical states of given quantum numbers J^*" 
and J: a one-particle state A (for example a or k) and a two-particle state P1P2 (for example 
TTTT or Kir). Let us study two eigenvalues of the 2x2 correlation matrix (for example using 
two interpolators from ([3])) 

Q,{t) = ±Tr[e-''^0,{t)0]m = ^ 5^(m|e-^(^-*)0,|n)(n|e-^*0j|m) (18) 



T 



and the relevant states n,m = A, P1P2, Pi, P2 giving [TT] 

Cij{t) = (0|C,|A)(A|C]|0)e-™'^* + (At|Oi|0)(0|Oj|At)e-™-4(^-*) (19) 

+ (P/|a,|P2)(P2|C]|P/)e-^^i(^-*)e-^^2* + (P2^|Ci|Pi)(Pi|C]|P2^)e-^^2(^-*)e-^^i* 

= {0\O,\A){0\Oj\A)* [e-^'At ^ ^-mMT-t)^ 

+ {0\Oi\PiP2){0\Oj\PiP2)* [e"^^i^2*±e^^^i^2(^-*)] 

+ (P||0,|P2)(P/|O,|P2)* [e-Ep^iT-t)^-Er,,t ^^-E,^iT-t)^-E,^t^ 

The signs "±" depend on the symmetry properties of the interpolators and states. We will 
assume "-I-" sign everywhere, since all our diagonal and non-diagonal correlators in the actual 
simulation are symmetric with respect to t ^ T — t (if the "— " sign occurs in some terms, 
the "backward propagation" described in Section Ilf of [27] may occur). 

If the coefficients in the third line were proportional to the coefficients in the second line 

(P/|Oi|P2) = Pi (0|Oi|PiP2) with Ri = R , independent of i , (20) 

the third line would depend linearly on the second line 

a,{t) = Y, zrzj*nt) (21) 

n=l,2 

with 

Zl = {0\O,\A), f\t) = e-"^* + e-'"^^^"*) , 

Z^ = (0|O,|PiP2), f\t) = e-^^1^2* + e-^^i^2(^-*) + R[e-Ep^iT-t)^-Ep,t ^ ^-.E^^iT-t)^-E^^ty 

(22) 

The eigenvalues of C{t)u^{t) = X^{t,to)C(to)u^(t) can be obtained in this case following the 
same steps as in [36] 

A"«) ^ ^ . (23) 

So, only if the condition (120|) would be exactly satisfied, the following applies: 
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Figure 12: Cosh- type effective masses ([T5|) for two eigenvalues of the generahzed eigenvalue 
problem with 2x2 correlation matrix (J19l20p and Ri ^ R2- A specific choice of the parameters is 
given in the text. 



(i) the eigenvalue corresponding to the one-particle state would have a time- dependence 
proportional to f^{t) and 

(ii) the eigenvalue corresponding to the two-particle state would have a time-dependence 
proportional to P{t). 

However, we believe that condition ( l20l) is not exactly satisfied, as argued after ( !T2|) in the 
main text. We have not rigorously derived the eigenvalues for this case. Based on the 
numerical example illustrated below, we however expect that eigenvalues have the form 



A"(t) = w" [e"^"* + e-^"(^-*)] + w" [e-'"^i*e-'^^2(^-*) + e-'"^2*e-'^^i(^-*)] 



(24) 



for two-particle states as well as for one-particle states. 

Let us demonstrate that on a specific example with tua = 0.6, mp^ = 0.3, rrip^ = 0.4, 
mp^P^ = 0.7, Zl = Zl = cos(0.2), Z^ = -Z\ = sin(0.2), T = 30, Ri = 2, i?2 = 0.4 and 
to = 4. We find the two eigenvalues for this case. Fig. [12] shows the corresponding cosh- 
type effective masses, defined in (IT5l) . The effective masses for the two-particle state P1P2 
as well as for the one particle state A drop near t ~ T/2. So none of them has the time- 
dependence e~^* + e~^^'^~*\ but a more complicated time dependence flMl) . The eigenvalue 
corresponding to the state A is proportional to e""^-^* + e"™^*-"^"*-* only in the special case 
Ri = R2- To demonstrate that the time-dependence of A"(t) is of the form ( l24l) to a very 
good approximation, we fit the resulting A"'(t) to ( 12^ with three free parameters En, w", w"'. 
The resulting En agree with input rriA or Ep-^p^ almost exactly and the resulting x^ of the 
filo is extremely small. 



11 



We attribute artificial error-bars ^"(t) = A"(i)/10 to the values of A"(t), which enter uncorrclatcd x^. 
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Appendix B: Discussion concerning omission of discon- 
nected diagrams in / = 0, 1/2 correlators 

In this appendix we show that the scattering state |/ = 2, /a = 0) cannot enter as an 
intermediate state in the connected correlation function {0^^^\0'''^^), while the scattering 
state |/ = |,/3 = I) cannot enter in (C^=V2|C)t^=V2). 

1=0 

The interpolators C[=°3 (^ 

Oi=% = -j^[2{dTu){uTd) + l{uTu){uTu) + l{dTd){dTd) - {uTu){dTd)] 

are |/ = 0, /a = 0) Clebsh-Gordon combinations of two 1 = 1 fields: |0, 0) = 4={2|1, 1)|1, — 1) — 
|1, 0)|1, 0)} = -^{-2{dTu){uTd) - ^[((irrf) - (uTu)]'^}, where we used fields qVq' with definite 
|/, J3) from Table [51 The | J = 2, Js = 0) state has flavor structure 

\2,0) = j^{2\lA)\l,-l) + 2\l,0)\l,0)} = j^{-2idT'u)iuT'd) + [idT'd)-iuT'u)]'}. 

We need to find out whether the state |2,0) can enter as an intermediate state in the 
connected part of the 1 = correlation function, therefore we evaluate the connected part 
of the matrix element 

(2, o\ol^=^)con = -^^ {2{dr'u){ur'd) - {ur'u){ur'u) - {dr'd){dr'd) + 2{ur'u){dr'd)\ (25) 

2{dTu){uTd) + l{uTu){uTu) + l{dTd){dTd) - {uTu){dTd)) 
= -^^ {AD{t) + [C{t) - D{t)] + [C{t) - D{t)] - 2D{t) + 2C{t) - 4C{t)} = 



Here D{t) and C{t) denote the "direct" and "crossed" connected contractions in the notation 
of [12], while six terms in the third line of (125 p refer to non-zero contractions of separate terms 
aa', bb', cc', dd', ad', da'; here abed refer to terms in the first row (I2S]) and a'b'c'd' refer to the 
second row. So we find that |2,0) state cannot couple to Ol^% via connected contractions. 

It is straightforward to show that {2,0\Ol^ )con = also for the diquark anti-diquark 
interpolators O^^^ ^. This is due to the cancellation between results from two terms in 
1 2, 0) oc {dT'u){uT'd) + {dV d){uT' u) + .. , which have the fiavor structure of Ol^^ ~ udud. 

So, the connected part of (2,0|(9^^^°) vanishes for all our interpolators and the |/ = 
2, /a = 0) state cannot enter as an intermediate state in our connected / = correlator. 

1=1/2 

Now let's see whether state |/ = |,/3 = \) 

\l \) = Tsilh -I)|l' 1) + v^li l)|l' 0)} = j=^{{-sT'd){dT'u) + {-sT'u)[{dT'd) - (uVu)]} 
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J 1 Icy 

can enter as an intermediate state of connected J = 1/2 correlator with 0]~ 
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0[~ll'^ = Yl (^rg)(grM) = {sTu){uTu) + {sTd){dTu) + (srs)(srM) . 

q=u,d,s 

The corresponding connected part of the matrix element is 

(i l\oY=l^')con = ^ {{sr'd){dr'u) - {sr'u){ur'u) + {sr'u){dr'd)\ 

{dTs){uTd) + {uTs){uTu) + {sTs){uTs)) 
= ^{D(t) + [C(t)-Z)(t)]-C(t)} = 

where three terms refer to non-zero contractions of separate terms aa', bb', ca'. So ||, ^) state 

cannot couple to Oi~ 3 via connected contractions. 

We also find that {^,^\Ol^ ) con = for diquark anti-diquark interpolators O^l (^. 
This is due to the cancellation between results from two terms in (|,|| oc {sT'd){dT'u) + 

{sT'u){dT'd) + .. , which have the flavor structure of 0^~^ ~ sddu. 

So, the connected part of (|, ^|0^'^^^^^) vanishes for all our interpolators and the |/ = 
|, I3 = |) state cannot enter as an intermediate state in our connected / = 1/2 correlator. 
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ensemble a [fin] ttIt^ [MeV] mT^a rriKCi N k, conf. 

C 0.1440(12) 318(5) 0.232(4) 0.391(3) 

B 0.1500(12) 469(4) 0.357(3) 0.462(3) 

A 0.1507(17) 526(7) 0.402(5) 0.465(3) 

Table 1: Two- flavor dynamical ensembles with Chirally Improved quarks [28]: pseudoscalar 
masses, Jacobi-smearing parameters {N, n) and the number of configurations are listed. All en- 
sembles have volume 16'^ x 32. 



777,^ [MeV] m-^a rrixa conf. 

230(7) 0.2332(56) 0.515(3) 300 

342(6) 0.3470(40) 0.545(2) 300 

478(8) 0.4840(32) 0.596(2) 300 

Table 2: The quenched simulation with overlap quarks is performed at two volumes 16^ x 28 and 
12^ X 28 at the same lattice spacing a = 0.200(3) fm |33j . The analysis is done at three pion masses 
above (these t?!,^ are determined at the larger volume). 
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I 


n 


m^ [MeV 


End 


w" 


w" 


to 


interp. 


tmin 


tmax 


xVdof 





1 


318 


0.36(4) 


0.33(8) 


-0.15(63) 


1 


C^12345 


8- 


-15 


0.015 





1 


469 


0.72(1) 


6.6(6) 


8.7(1.0) 


3 


C^12345 


8- 


-15 


0.0065 





1 


526 


0.81(1) 


8.7(7) 


9.2(9) 


3 


C^12345 


8- 


-15 


0.14 





2 


318 


0.41(2) 


0.061(10) 


0.049(16) 


1 


Cl2345 


8- 


-15 


0.088 





2 


469 


0.72(1) 


3.8(3) 


5.3(4) 


3 


Cl2345 


8- 


-15 


0.024 





2 


526 


0.82(1) 


5.3(5) 


6.6(5) 


3 


Cl2345 


8- 


-15 


0.018 





3 


318 


0.93(9) 


0.23(12) 




1 


C^12345 


6- 


-9 


0.36 





3 


469 


1.11(5) 


10.5(3.7) 




3 


C^12345 


7- 


-9 


0.28 





3 


526 


1.07(6) 


7.7(2.9) 




3 


C^12345 


6- 


-10 


1.1 


2 


1 


318 


0.45(2) 


0.33(4) 


0.32(8) 


1 


Ol23 


8- 


-15 


0.0086 


2 


1 


469 


0.736(5) 


0.47(2) 


0.54(3) 


1 


£^123 


8- 


-15 


0.016 


2 


1 


526 


0.827(7) 


0.57(4) 


0.54(4) 


1 


£^123 


8- 


-15 


0.0084 


2 


2 


318 


1.33(19) 


0.73(91) 




1 


(^123 


6- 


-9 


0.0039 


2 


2 


469 


1.33(4) 


0.41(13) 




1 


C^123 


7- 


-10 


0.25 


2 


2 


526 


1.35(8) 


0.40(27) 




1 


C^123 


8- 


- 11 


0.60 


1/2 


1 


318 


0.55(3) 


2.8(6) 


1.5(1.4) 


3 


Cl2345 


8- 


-15 


0.077 


1/2 


1 


469 


0.826(8) 


8.1(6) 


9.2(9) 


3 


Cl2345 


8- 


-15 


0.032 


1/2 


1 


526 


0.87(1) 


9.4(9) 


9.4(9) 


3 


Cl2345 


8- 


-15 


0.11 


1/2 


2 


318 


0.636(8) 


3.4(2) 


3.2(3) 


3 


C^12345 


8- 


-15 


0.011 


1/2 


2 


469 


0.852(5) 


6.7(3) 


7.2(4) 


3 


C^12345 


8- 


-15 


0.048 


1/2 


2 


526 


0.907(7) 


8.3(6) 


7.4(4) 


3 


C^12345 


8- 


-15 


0.012 


1/2 


3 


318 


1.05(7) 


8.8(40) 




3 


Cl2345 


7- 


-10 


0.11 


1/2 


3 


469 


1.16(6) 


9.9(45) 




3 


Cl2345 


8- 


-10 


0.57 


1/2 


3 


526 


1.14(4) 


7.9(23) 




3 


C^12345 


7- 


-10 


0.32 


3/2 


1 


318 


0.630(9) 


0.36(2) 


0.34(4) 




0123 


8- 


-15 


0.0088 


3/2 


1 


469 


0.844(5) 


0.50(2) 


0.55(3) 




C^123 


8- 


-15 


0.042 


3/2 


1 


526 


0.896(6) 


0.60(4) 


0.57(4) 




C^123 


8- 


-15 


0.013 


3/2 


2 


318 


1.27(9) 


0.33(21) 






C^123 


7- 


-10 


0.62 


3/2 


2 


469 


1.40(5) 


0.50(19) 






C^123 


8- 


-10 


0.017 


3/2 


2 


526 


1.34(9) 


2.5(2) 






C^123 


9- 


-11 


0.058 



Table 3: Extracted energies End (a^^ ~ 1.3 GeV) together with w"^ in the dynamical simulation 
for all isospins. The fit form ()13p is used whenever Wn is provided, while fit form (J17p is used where 
w'^ is not provided. The interpolator basis, to, fit ranges and uncorrelated x^/(degrees of freedom) 
are also presented. The n = 1 states with / = 0, 1/2 and the lowest ttItt have badly determined tD" 
since they have almost flat cosh-type effective mass, which indicates they are roughly consistent 
with (llTp and iD ~ (finite T effect is less significant at low m.^ |llj). 
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I 


n 


m^ [MeV 


E^a 


w" 


w" 


to 


interp. 


<^m.in 


i^max 


xVdof 





1 


230 


0.40(3) 


1.2(6) 




3 


Ol23 


8- 


-10 


0.0013 





1 


342 


0.72(2) 


3.1(6) 


1.5(4) 


3 


C^12345 


8- 


-13 


0.0010 





1 


478 


1.03(2) 


7.2(14) 


2.0(3) 


3 


Cl245 


8- 


-13 


0.0020 





2 


230 


0.45(4) 


0.15(6) 




3 


Ol23 


8- 


-10 


0.039 





2 


342 


0.75(8) 


0.45(44) 


0.57(9) 


3 


C^12345 


8- 


-13 


0.0084 





2 


478 


1.03(6) 


0.46(43) 


0.38(4) 


3 


C^1245 


8- 


-13 


0.018 





3 


230 


1.05(21) 


0.8(12) 




3 


Ol23 


6- 


-8 


0.41 





3 


342 


1.15(6) 


1.3(5) 




3 


C^12345 


7- 


-9 


0.018 





3 


478 


1.33(5) 


1.2(5) 




3 


C^1245 


8- 


-10 


0.053 


2 


1 


230 


0.54(9) 


3.3(19) 


2.5(6) 


3 


C^123 


8- 


-13 


0.085 


2 


1 


342 


0.719(9) 


1.8(2) 


1.64(9) 


3 


C^123 


8- 


-13 


0.0027 


2 


1 


478 


1.032(7) 


2.8(2) 


1.90(7) 


3 


Ol23 


8- 


-13 


0.16 


2 


2 


342 


1.00(37) 


0.05(55) 




3 


Ol23 


7- 


- 10 


0.21 


2 


2 


478 


1.50(9) 


0.62(52) 




3 


Ol23 


8- 


-11 


0.21 


1/2 


1 


230 


0.76(3) 


2.7(10) 


1.6(5) 


3 


C^12345 


9- 


-13 


0.00034 


1/2 


1 


342 


0.94(2) 


4.5(8) 


1.7(3) 


3 


Cl2345 


9- 


-13 


0.0030 


1/2 


1 


478 


1.16(2) 


8.1(13) 


1.9(2) 


3 


C^12345 


9- 


-13 


0.0038 


1/2 


2 


230 


0.76(3) 


1.1(3) 


1.1(2) 


3 


C^12345 


9- 


-13 


0.0030 


1/2 


2 


342 


0.94(2) 


1.2(2) 


0.90(6) 


3 


C^12345 


9- 


-13 


0.00018 


1/2 


2 


478 


1.18(1) 


2.0(2) 


0.93(4) 


3 


C^12345 


9- 


-13 


0.029 


1/2 


3 


230 


1.23(7) 


1.7(8) 




3 


C^12345 


7- 


-10 


0.014 


1/2 


3 


342 


1.29(9) 


1.0(7) 




3 


C^12345 


8- 


-10 


0.096 


1/2 


2 


478 


1.47(4) 


1.7(5) 




3 


C^12345 


8- 


-10 


0.22 


3/2 


1 


230 


0.771(9) 


2.1(2) 


2.2(1) 


3 


£^123 


9- 


-13 


0.0062 


3/2 


1 


342 


0.937(6) 


2.3(1) 


1.93(7) 


3 


£^123 


9- 


-13 


0.0047 


3/2 


1 


478 


1.156(5) 


3.45(2) 


2.15(7) 


3 


C^123 


9- 


-13 


0.046 


3/2 


2 


230 


1.19(29) 


0.09(22) 




3 


C^123 


8- 


-10 


0.00014 


3/2 


2 


342 


1.27(11) 


0.09(12) 




3 


C^123 


9- 


-12 


0.23 


3/2 


2 


478 


1.55(7) 


0.47(33) 




3 


Ol23 


9- 


-11 


0.12 



Table 4: Analogous to Table El but for the quenched simulation at the larger volume V = 16^ x 28 
and a~^ ~ 1 GeV. 



/,/3) 


1,1) 


1,-1) 


1,0) 


- -) 

2' 2/ 


i _i\ 

2' 2/ 


qTq' 


dVu 


—uTd 


-j^[dTd-uTu 


sTu 


sTd 



Table 5: qTq' with definite |/, /s), which are used to build two-particle interpolators via Clebsh- 
Gordan coefficients in Appendix B. 
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